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We study atomic autoionization processes in the time domain. With the emerging attosecond 
extreme vacuum ultraviolet and soft x-ray pulses, we first address how to characterize the time 
evolution of the decay of a discrete state into a degenerate continuum. A short pump beam generates 
a number of resonance states in a series and the nearby background continuum, and the resultant 
wave packet evolves with time until the full decay of the bound states. Taking the 2pns( 1 P°) 
resonance series embedded in the 2sep( 1 P°) continuum in beryllium atom as an example, the time 
evolution of the autoionizing wave packet in energy domain and in coordinate space is calculated and 
analyzed, where Fano profiles build up in the photoelectron energy during the process. A proposed 
pump-probe scheme assumes that the probe beam ionizes the 2s inner electron in the wave packet. 
The lifetimes of the resonances and the photoelectron energy distribution can be obtained from the 
ionization yield versus the time delay of the probe. 

PACS numbers: 32.80.Zb,32.80.Fb 



I. INTRODUCTION 



Scattering processes - photoionization, electron colli- 
sion, etc. - have been important tools for studying phys- 
ical systems since the early days of the development of 
atomic theories. Autoionization, first observed by Bcut- 
ler [l[ in the photoabsorption of rare gas atoms, is a 
phenomenon where the energy spectra is characterized 
by asymmetric resonance peaks, which were explained 
by Fano in his seminal paper in 1961 0. In this pa- 
per, asymmetric resonance in photoabsorption occurs in 
the energy region where a discrete state is embedded in 
the continuum. Photoabsorption populates both the dis- 
crete state and the continuum states. The interaction 
of the discrete state with the continuum results in the 
"autoionization" of the former, where the discrete state 
interferes with the directly populated continuum states 
to form asymmetric Beutler-Fano profile. Today Fano 
resonances have been used to describe a wealth of physi- 
cal phenomena in atomic, molecular and optical systems, 
condensed matter systems 0,0], nanostructures[|[ as well 
as many other fields. Essential to Fano's theory is that 
each resonance can be characterized by three parameters, 
the "unshifted" resonance position E r , the shape param- 
eter q and the resonance width T. The resonance width 
characterizes the decay time or the lifetime of the discrete 
state. Since the typical resonance width of an atomic or 
molecular system is in the order of 10 _1 eV, which corre- 
sponds to the lifetime of as short as a few femtoseconds, 
almost all measurements on these resonances have been 
carried out in the energy domain, using high-resolution 
spectroscopy. Today a large set of experimental and the- 
oretical data of these resonances have been collected, and 
they provide important structural information of the sys- 
tems. 

In the past decade, with the advent of laser technology, 
single attosecond pulses (SAP) are becoming available [6|. 
Such pulses opened up a new research platform for study- 



ing electronic dynamics, such as the autoionization dy- 
namics. Indeed, the first SAP experiment in atomic sys- 
tem was to determine the lifetime of the Auger decay of 
an M -shell hole in Kr created by extreme vacuum ul- 
traviolet (XUV) attosecond pulse by Drescher et al 0]. 
For Auger decay, the resonance profile is mostly sym- 
metric, and thus only the decay width or the lifetime of 
the inner-shell hole was determined. In Ref. Q, the life- 
time is determined by creating the inner-shell hole in the- 
presence of an IR laser. By changing the time delay be- 
tween the XUV pulse with the IR, photoelectron spectra 
can be used to "deduce" the lifetime of the inner-shell 
hole - based on the approximate "streaking" theory Q. 
Such streaking theory has been generalized to Fano res- 
onances by Zhao and Lin as well as by others [13, [H| ■ 
However, the photoelectron spectra from such XUV+IR 
setup are very complicated and only the lifetime of the 
Fano resonance can be extracted. 

Due to the insufficient fiuence of the SAP's so far, ex- 
isting measurements using attosecond pulses for Auger 
decay or autoionization have all been carried out using 
the XUV+IR pump-probe setup where the two pulses ac- 
tually overlap in time@,[ll,[ill- Rigorously speaking, the 
IR is not probing the resonances created by the attosec- 
ond XUV pulse in these experiments since the photoelec- 
tron is created in the IR field. With improved technology 
in attosecond pulse generation and the recent advances 
in X-ray free-electron lasers (XFEL), we envision that in- 
tense attosecond pulses will become available in the near 
future, such that a true time-dependent investigation of 
the autoionization process will be possible. In this arti- 
cle, we analyze the time evolution a Fano resonance, or 
the autoionization dynamics of the system in the field- 
free condition before the discrete state fully decays into 
the continuum. For example, if the discrete state has a 
lifetime of 20 fs while the pump pulse is only 1.5 fs, we 
ask what is the time evolution of the wave packet asso- 
ciated with the decay process after the pump pulse, and 
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how the energy profile eventually evolves into what is 
characterized by the Fano resonance shape. 

We comment that the time evolution of a Fano reso- 
nance has been studied in several calculations by solving 
the time-dependent Schrodinger Equation of helium or a 
model atom H^-lljj. with or without the IR probe. 
These studies did not address how to retrieve the time 
evolution of the resonance profiles. Similarly, autoioniza- 
tion of double Rydberg states in the time scale of picosec- 
onds has been studied [171] . In the meanwhile, the interac- 
tion of a doubly excited state with a whole Rydberg series 
also falls into the similar category [lj|[l9[ . Understanding 
the time evolution of the two-electron wave packets and 
the issue of how to probe their evolution lie at the heart 
of understanding and controlling two-electron dynamics. 

In this article, we use two pulses that are separate in 
time. The probe is always placed strictly after the pump 
so that the field-free evolution between the pulses can be 
defined. The physical process is understood roughly as 
the decaying of the bound state to the continuum due to 
the interaction between them. In Sec. HH starting with 
Fano theory which is for stationary states, we construct 
the time-dependent wave function of the system which 
describes the time evolution of the electron wave packet 
after the pump pulse. The wave function is recasted in 
the basis set of the bound state and the continuum states, 
where it contains the time information for the bound and 
the continuum components. The wave packet for the 
continuum part is shown to evolve to the typical Fano 
resonance profile at long time. We thus define the con- 
tinuum wave packet at short time as the time-dependent 
Fano profile of the autoionization. Since the pump pulse 
is short, it inevitably spans a broad-band in the energy 
spectra. Assuming that the background energy spectrum 
generated by the pump is described by a Gaussian dis- 
tribution, we generalize the theory to include all the res- 
onances within the bandwidth. In Sec. IIIII we apply our 
theory to the resonances in beryllium to investigate wave 
packet dynamics after the pump pulse. In Sec. IIVI we 
demonstrate how to probe the time-dependent Fano res- 
onances, using the example in beryllium given in Sec. IIIII 
A short summary and future outlook is given in Sec. |Vj 



II. THEORY 



A. Fano's theory of resonance 



In a multi-electron atomic system, the high-lying dis- 
crete states are embedded in the continuous spectrum if 
these discrete energies are higher than the binding energy. 
Both the bound and continuum states are described in 
terms of configurations. According to Ref. Q, with one 
bound state and its nearby continuum, the system is gov- 



erned by the total Hamiltonian given by 

(a\H\a) = E r , 
(j3 E \H\a)=V E , 
{(3e>\H\(3 e ) = ES(E' - E), 



(1) 

(2) 
(3) 



where |a) and \@e) are bound and continuum configura- 
tions respectively. The off-diagonal terms Ve represent 
the mixing strength between the bound and the contin- 
uum. The diagonal terms E r and E are bound and con- 
tinuum energies respectively. 

The eigenstates \i/je) of the atom are obtained by di- 
agonalizing the Hamiltonian in Eq. QJ3] and are given by 
Fano as 
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Here we assume that V E = V is a real constant in the 
narrow energy region for each Fano resonance, provided 
that the bound and the continuum states are in real rep- 
resentation. 

Originally Fano's theory was applied to photoabsorp- 
tion processes. The transition amplitude is ce = 
(ipE\T\g) where \g) is the ground state, \tp E ) is the excited 
state of energy E, and T is the photon transition opera- 
tor. The states \g) and \4>e) are eigenstates of the atomic 
Hamiltonian H 7 so the coefficients c E do not change with 
time. In Ref. Q, the continuum wave functions are taken 
to be real standing waves. The absorption cross section 
is proportional to |ce| 2 , which contains the features of 
Fano profiles. Until recently, all photoabsorption mea- 
surements are performed using light sources of hundreds 
of picoseconds or longer, where they can be considered 
as monochromatic. Thus, \iPe) of a single E is pro- 
moted in each measurement. To trace out a Fano res- 
onance, monochromatic lights are tuned across the res- 
onance while the cross section for each energy point is 
recorded. 

Although Fano's theory was applied to photoabsorp- 
tion originally, it can be trivially generalized to study 
photoelectron angular distributions. While the contin- 
uum functions in Eq. 3] arc taken to be real standing 
waves, \iI>e) is an eigcnstatc of the Hamiltonian. To ob- 
tain the scattering amplitude for an electron with mo- 
mentum fc, one only needs to project the standing wave 
onto the momentum eigenstate k of the photoelectrons, 
in energy-normalized form given by 




^ Ul {kr)Y lm {f)YL{k). (8) 
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Here ui (kr) are the radial wave function of the electron 
taken as real standing waves as in Fano's theory, r\i are 
the total scattering phase shifts, and the Y's arc spherical 
harmonic functions. 

Following the notation in Ref . , the width of the res- 
onance is defined by T = 2irV 2 . Assuming that (/3e|T|(?) 
is energy- independent near the resonance energy E r , the 
q-paramctcr is defined by 



(a\T\g) 



nV(p E \T\g) 



(9) 



Without the flat-background and the constant V assump- 
tions, | a) will be replaced by the "modified bound state" 
defined in Ref. Q which contains a term based on un- 
even Ve and \f5 E ) near E r . We assume this modification 
is negligible in general cases and consider the g-parameter 
as defined in Eq. [HI Fano's theory predicts the resonance 
profile in a simple form, where the associated parameters 
E r , r, and q can be extracted by fitting the form to the 
measured profile. 

With the advent of available light pulses of a few fem- 
toseconds or shorter, one can use such a short pulse to 
ionize the atom. After the pump pulse is over, the time- 
dependent wave function of the atom is given by 
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where in a single pulse, the integral in Eq. [10] covers a 
bandwidth more than the width of a resonance. Such a 
resonance can be measured, long after the pump, in a 
high-resolution spectrometer. This amounts to project- 
ing out the time-dependent wave function \^(t)) over a 
long time. Thus the whole Fano resonance is generated 
in a single pulse and its profile is given by \ce\ ■ 

We now consider the situation where the pump pulse is 
shorter than the lifetime of a Fano resonance. One would 
like to "see" the dynamics of the autoionization. How- 
ever, such information is not contained in \ce\ 2 since \iPe) 
are eigenstates of the atomic Hamiltonian and the coef- 
ficients ce are constant when the pump field is absent. 
Thus in Eq. 1101 the time information is hidden and not 
recognized, even though we "know" that autoionization 
is happening! 



B. Time-dependent wave function in configuration 

basis 



It is convenient to define the excited wave function 
\ty ex (t)) as the total wave function \^(t)) excluding the 
ground state part since the ground state does not partic- 
ipate in autoionization. Instead of using eigenstates as 
expressed in Eq. [TU1 we alternatively expand the excited 
wave function in configuration basis as 

|*ex(*)> = d a (t)\a) + ( d E (t)\p E )dE. (11) 



Since |a) and \p E ) are not eigenstates of the atomic 
Hamiltonian, the coefficients d a {t) and dsit) are time- 
dependent. One expects that at long time (compared to 
the lifetime), the coefficient d a (t) will go to zero and all 
the information of the total wave function will be con- 
tained in dE{t). However, the time evolution of these co- 
efficients cannot be measured in laboratory - they change 
in the time scale of femtoseconds or less, while laboratory 
electronics can only measure changes of the order of pi- 
coseconds or longer. Here we define \dE(t)\ 2 as the "time- 
dependent Fano resonance profile", or in short, "reso- 
nance profile". Below we discuss how these coefficients 
behave. The measurement issues will be addressed later. 

Assuming that the initial values of d a (t) and dE{t) 
are known, with the aid of Eq. |4][7l the time-dependent 
coefficients are solved and given by 
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where the functions g E (t) and fEE'it) are defined by 
V 



g E (t) = 

fEE'it) = 
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The superscript (0) denotes the initial values, i.e., 

4 0) = («|*(0)), 
4 0) = (/^l*(0)). 



(14) 
(15) 



(16) 
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The solutions in Eq. [T2lfT5l are calculated exactly for a 

i E , E r , and T, where 



given set of parameters da\ d^ 



these parameters can be either obtained by running a 
separate program such as Time-dependent Schrddingcr 
Equation (TDSE) code or time-dependent perturbation 
calculation, based on the prior knowledge of the pump, 
or extracted from the experimental values in the time- 
integrated spectrum. In the latter case, with q defined 
by Eq. [9] and rewritten as 



d {0) 



irVd 



CO) 



(18) 



E=E,, 



the initial coefficient d^ can be extracted since the value 
of q and the background magnitude of can be ob- 
tained from the Fano profile after long time (the phase 
of djp can be deduced if photoelectron angular distribu- 
tions are measured) . Equations [T2"l[T5l are the backbone 
of our model. 

Note that the initial coefficients d^ and djp are in 
general complex values while q is real. Within the narrow 
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range of a resonance width, the energy E in dg is close 
to the bound state energy E r such that the phases of da^ 
and d!^) are assumed to be the same in the neighborhood 
of the resonance, and the corresponding q, representing 

the ratio between da and a§ , is real. When d^ spans 
a wide energy range covering many resonances, its phase 
varies smoothly with E, but in each resonance region, 
the statement above is still true, and each q is real. 

The decay lifetime of the bound state is defined as 
T = 1/r. Rigorously speaking, the dynamics of the au- 
toionization is not only about the decay of the bound 
state, but also its coupling to the continuum (See Eq . fTTT) . 
The wave packet can propagate from the continuum to 
the bound or the other way around, or between contin- 
uum states of different energies through the bound state. 
Thus, the lifetime T, determined by the interaction ma- 
trix V, characterizes the time scale of the whole event. 

We can generalize Eq. fTWl5l to include multiple res- 
onances by applying Fano's theory for many bound 
states and assuming that higher-order perturbation 
by V is negligible. The excited wave function is 

* e *(i) = X) d "(*)l a ">+ f d E (t)\0 E )dE, (19) 

n J 

with the coefficients given by 
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(20) 



(21) 



where the label n is the index for the nth bound state. 
The ^-function and the /-function, now subscripted by 
n, are calculated independently for each resonance with 
the associated parameters as given by Eq. [1411151 

To simulate the effect of the finite range of V while 
keeping the calculation simple, the profiles can be fur- 
ther adjusted by multiplying each "resonance term" in 
Eq. [5T] (in the {} bracket) by a Gaussian window func- 
tion to limit the range of each resonance. For an isolated 
resonance, the effect of the window function is minimum 
or unnecessary because the window function mainly re- 
moves "the wings" but not the center of the resonance, 
where the central part is the significant and informative 
part. On the other hand, when dealing with many res- 
onances, without the window functions, the wings may 
be overly extended and unnaturally perturb the nearby 
resonances. In the latter case, multiplying the window 
functions is an effective and efficient solution. 



d^ = const. This initial condition, dubbed ICi, is valid 
near the resonance positions if the pump bandwidth is 
broad enough relative to the resonance widths so that 
the populated profile is nearly fiat around E r . For con- 
ceptual presentation, only the single-resonance formulae 
Eq. [T2TT51 are referred here, while the idea is the same 
for the many-resonance case. Defining the scaled time s 
by s = Ft and scaled energy e by e = 2(E — E r )/T, the 
bound and the continuum coefficients under ICi are 



d a ( S ) = <4°)(l--)e-h 
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d e (s) = d^- 1 
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(q 



-(q-i)e ^ 



where in the scaled energy, q-parameter becomes 
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(22) 
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and the factor e~' lErt disappears because e is measured 
relative to the resonance energy. While the coefficients 
given by Eq . [2211231 represent the evolution for s > 0, their 
limits at s — > arc different from the initial coefficients 
da^ and di°^ that we have defined at s = 0. This discon- 
tinuity means an immediate jump of the wave function 
right after the evolution starts. This unphysical behav- 
ior is due to the integral of g-function in Eq. [TJ] and /- 
function in Eq. ll5l over infinite energy range at very small 
times, where the contributions from E much higher or 
lower than E r are misleadingly counted in, while in prac- 
tice, only the energy range comparable to T should be 
considered. Nevertheless, the large-time behavior is reli- 
able since both g- and /-functions are compressed more 
and more into the near-resonance region with increasing 
time. 

In the s — > 00 limit in Eq. [2211231 the bound state 
decays to zero, but the continuum recovers the general 
form of Fano profile, 



lim K( s )| 2 = |d(°)| 
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(25) 



We emphasize that Eq. [25] is a mathematical conse- 
quence. It cannot be obtained from the measurement 
since continuum configurations are not eigenstates of the 
total atomic Hamiltonian. In the special case where 
d^ = 1 and di°^ = 0, the system starts with only bound 
state, and Lorentz profile will be recovered as s — > 00, 



i.e. 



lim \d e {s)\ 2 = - 



1 1 



1 



(26) 



1. Flat initial continuum (IC) distributions 

To study the qualitative behavior of the wave func- 
tion, we first assume that the initial profile is flat, i.e. 



Fig. [T] shows the spectra at different times for q — > 00, 
q = 2, q = 1, and q = 0. Each sub- figure demonstrates 
how its spectrum morphs into Fano profile. A feature 
found in the profile is the oscillation in the wings during 
the evolution, which shrinks its energy period when time 
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The initial coefBcicnts are given by 




FIG. 1: The resonance profile |d £ (s)| 2 calculated by Eq. l23l for 
(a) q ->■ oo, (b) q = 2, (c) q = 1, and (d) q = 0. For (b)-(d), 
\dl | 2 is normalized to 1. In all cases, the profile builds up 
with time with small oscillations in energy before reaching the 
final Fano profile. 



dn = IQnKVn exp 



1 f E n - uj 
2 



D 



d E = 7 exp 



1 (E-oj 

2 \ D 



where 



= y/nD + iln-n-Vn) 2 exp 



E n - 
D 



(27) 
(28) 

(29) 



This initial condition is dubbed IC2- The initial contin- 
uum state has a Gaussian energy distribution centered at 
u> with standard deviation D (FWHM = 1.665D). The 
initial coefficients are determined by the pump width and 
the resonance parameters E n , T n (or equivalently, V n ), 
and q n . For D >• T n , IC2 is reduced back to ICi near the 
nth resonance. The comparison between the profiles un- 
der the ICi and IC2 is plotted in Fig. [2] for q=l, where 
the bandwidth (FWHM) for IC 2 is 15T. As seen in Fig. 03 
with the modification from ICi to IC2, the s = profile 
changes from a flat background to a Gaussian shape, and 
the s = 1 and s = 2 profiles change similarly. However, 
the featuring structures are preserved. The modification 
should be slight if the bandwidth is much wider than F, 
which is true for a sufficiently short pump. 



increases, and disappears as s — > 00. This oscillation, due 
to the phase between the two exponential terms in Eq.l23[ 
is not seen in the conventional time-integrated spectra. 
Although analytically simple, the ICi simplification is 
conceptually rich and offer an easy visualization for the 
autoionization process. 




2. Gaussian initial continuum distributions 



The condition ICi is valid only in the vicinity of a 
resonance and not applicable to the whole energy range. 
When a short pump pulse is used to initiate the system, 
one should cover its whole bandwidth and the state vec- 
tors must be properly normalized. Here we consider the 
normalization where the total probability for the bound 
component \d n (t)\ 2 plus the continuum component 
J \d,E(t)\ 2 dE is set to 1.0, where the ground state \g) is 
not considered since it is not involved in the autoioniza- 
tion. This normalization is automatically fulfilled by the 
exact solution once the initial coefficients are normalized. 
Additionally, we assume that the pump pulse has a Gaus- 
sian envelope, resulting in a Gaussian profile in energy. 



FIG. 2: The resonance profile for ICi and IC2 initial con- 
ditions. The solid curves are carried out with the IC2 initial 
condition defined by q = 1 and the FWHM bandwidth of 15F, 
or D = 9r. The dashed curves are carried out under ICi. 



3. Growth and decay of a bound state 

Considering an isolated resonance, in ICi, the total 
bound state is an exponential decay function of time as 
described by Eq. 1221 However, under the normalization 
condition of IC2, which is more realistic, the bound state 
decay is not necessarily monotonic. For the IC2 solu- 
tion, the bound state coefficient in Eq. [T^] and its two 
terms are shown in Fig. [3] in their absolute squares. The 
two panels arc for two different initial bandwidths of the 
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continuum centered at the resonance. The first term, 
which comes from the initial bound state, is always an 
exponential decay. The second term, which comes from 
the initial continuum, goes through a rising period up to 
about s = 0.2 in Fig.^a) and s = 0.1 in Fig.^b) before 
it decays. This rising period will shorten if the initial 
bandwidth increases or the pump duration shortens. In 
the limit of infinite bandwidth, both the first and the 
second terms will be monotonic decay functions, and the 
solution in ICi will be recovered. Note that even in a sit- 
uation where the length of the rising period is comparable 
to the lifetime, such as what Fig. G^a) shows, the total 
bound state population will eventually be dominated by 
the exponential decay. This growth part of the bound 
state is usually much shorter than the decay lifetime, so 
to "watch" this detail structure, a higher temporal reso- 
lution is required. 
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FIG. 3: The bound state coefficient defined by Eg. 1 121 and its 
compositions for q = 1 and (a) D — 10T, and (b) D = 20T, 
where D is the standard deviation of the pump beam (See 
Sec lII B 3| . When D increases relative to T, the rising portion 
of the second term at the beginning of time gets shorter, and 
in the limit J)»r, Eg. [221 (ICi ) is recovered. 

In the special case for d!^ = 0, the pump pulse does 
not directly populate the bound state of the system. The 
q-parameter defined by Eq. [T5] is 0. However, this defi- 
nition is based on the assumption that V and the con- 
tinuum background are constant near the resonance, and 
it excludes the modification term for the bound state. 
In an actual situation where no original bound state \a) 
was populated, the theory needs to be modified to in- 
clude higher order terms. Nonetheless, even if we neglect 
the initial population of the original and the modified 
bound state, the time-dependent probability as given by 
Eq. [12] will still have contribution from the initial con- 
tinuum population . In other words, we will still 
have the "decay" of the bound state even if it was not 
populated by the pump pulse. This somewhat awkward 
result is best understood in the following way. First re- 
call that the initial bound state population depends on 
the transition operator connecting it to the ground state, 
while the "configuration interaction" V which connects 
the bound state to the continuum is always present. This 
"configuration interaction" feeds the bound state from 
the continuum, but the continuum states also draw the 



bound state component and take it to large distances. 
If these continuum states were substituted by highly ex- 
cited bound states, the outer electron would return to 
recollide with the inner electron(s) and repopulate the 
bound state. This would correspond to situations stud- 
ied in Ref. (l7l - fl9| . To initiate the oscillation, the bound 
state does not have to be excited by the pump pulse. The 
ever present " configuration interaction" is the one which 
is responsible for the oscillation. Such "interaction" is 
the property of the Hilbert space but not of the specific 
excitation mechanism. 



III. APPLICATION TO RESONANCES IN BE 

In atomic systems, autoionization is understood only 
in the multiple-electron picture where electron-electron 
correlation is considered. Our theory is applied here to 
calculate the autoionizing wave packet of the 2p4s( 1 P°) 
resonance and of the 2pns( 1 P°) (n = 3 to 9) resonance 
series embedded in the 2sEp( 1 P°) continuum in beryl- 
lium atom, generated by the photoionization from the 
2s 2 ( 1 S e ) ground state. For a short pump, we assume 
that the initial continuum distribution can be approxi- 
mated by a Gaussian distribution. Only the two outer 
electrons in beryllium are assumed active, and the Is 2 
core electrons are frozen. The bound orbitals are Slater 
type orbitals where their analytical forms are determined 
by fitting the numerical Hartree-Fock calculations. The 
continuum basis functions are calculated with the model 
potential given in the literature [20j for I = 1 stand- 
ing waves 2l[. The resonance parameters are taken 
from the experiment of Wehlitz et al [22j], where E r is 
2.789 eV above the ionization threshold, q = —0.52, and 
T = 0.174 eV, or T = 3.78 fs. The pump beam has the 



carrier frequency right at the resonance and the duration 
of 2 fs, or mean bandwidth of 0.912 eV. 

For the 2pAs resonance case, the excited wave function, 
Eq. 111! is given by 

+ I d E ,(t)q> 2sE , p (ri,r 2 )dE', (30) 



where the (^-functions are constructed by corresponding 
bound orbitals and continuum waves, and symmetrized 
between the two electrons. The one-electron density is 
defined by 



P(ri,t) 



^(rl,r1;i)| 2 rfr 2 dr!idr! 2 dr 2 , (31) 



where the integral sums over the angular dependence as 
well as the radial part of one of the electrons. In the large- 
r region, p(r±,t) represents the electron density of the 
autoionizing electron. In this section, energy, time, and 
distance are in electron Volts (eV), femtoseconds (fs), and 
Bohr radius (a.u.) respectively unless otherwise specified. 
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The bandwidth in energy and the duration in time are 
defined in FWHM. 

First we show how the wave packet moves in the coor- 
dinate space. Fig. 2(a) displays the electron density at 
t = 0, 2 fs, and 4 fs, where the time scale is about the 
lifetime of the resonance T = 3.78 fs. Before the propaga- 
tion starts (t — 0), a dominant spatial distribution covers 
0-60 a.u. This distribution is mainly contributed by the 
continuum, where the percentage of initial bound state is 
as low as 7%. The spatial oscillation with period of about 
10-15 a.u., or wave number of about 0.4-0.6 a.u., matches 
the energy distribution centered at the resonance energy 
of 2.789 eV. As time goes on, the electron wave packet 
moves outward toward large r. When time is long com- 
pared with the autoionization lifetime, the motion of the 
electron wave packet can be analyzed more easily. In 
Fig. [4(b), the long-time behavior of the wave packet are 
plotted for t = 40 fs, 80 fs, 120 fs, and 160 fs. The 
electron is always confined in some spatial range, but it 
gradually spreads wider when time increases. By recog- 
nizing and tracing the dominant two peaks in Fig. [4(b) 
over time, we find that the corresponding velocities are 
almost constant. It is similar to a free electron wave 
packet traveling in space, where its average velocity keeps 
the same, but its spatial distribution spreads. In this long 
time regime, the electron moves in a very predictable way 
to the large distance toward the detector. 

Next, we discuss the evolution of the resonance pro- 
file \d,E(t)\ 2 , i.e., the time-dependent Fano profile, for the 
same wave packet in Eq. 1301 In Fig. [5(a), the photoelec- 
tron energy profile is shown at t = 0, 1.5 fs, 3 fs, and 
4.5 fs, i.e., near the lifetime of the resonance T = 3.78 fs. 
The profile starts evolving with a Gaussian shape at 
t = determined by the pump pulse, and changes quickly 
within its lifetime. At t = 4.5, just a little longer than 
the lifetime, the double-peaked shape characteristic of 
autoionization can be seen, but the shape is not quite 
"settled" into the Fano profile. The long-time behavior 
is shown in Fig. [5(b). At t = 10 fs to 20 fs, the resonance 
profile is already quite close to the final Fano profile that 
can be determined in the standard energy domain mea- 
surements. It is also physically sensible that the electron 
velocities observed in Fig. [4(b) can be mapped onto, al- 
though not exactly, the energies of the two highest peaks 
in Fig. 15(b). 

We have discussed an isolated resonance. Practically, 
a short pump pulse is likely to populate more resonances 
at once. Here we consider the case including the 2pns 
resonances up to 2p9s in Be. Even though the 2p2s state 
(about 5.277 eV above the ground state (23|) is lower 
than the 2s threshold and will not show up in the con- 
tinuum profile, it can in general interact with the con- 
tinuum beyond the 2s threshold and is included in the 
present calculation. The parameters for these resonances 
are generated by setting [i = 0.6, T x v z = 0.71 eV, 
and q = —0.8, where [i is the quantum defect and v 
is the effective quantum number. The first few widths, 
starting from 2p3s, are 0.514, 0.181, 0.0833, and 0.0451, 
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FIG. 4: The electron density for the 2p4s resonance in Be by 
a pump pulse of 2 fs and central frequency at the resonance. 
The decay lifetime of the resonance is 3.78 fs. The electron 
density is shown for (a) short-time behavior at t = 0, 2 fs, 
and 4 fs, and for (b) long-time behavior at t = 40 fs, 80 fs, 
120 fs, and 160 fs. 



in eV's, and the decay lifetimes are 1.28, 3.64, 7.90, 
and 14.6 in femtoseconds, respectively. The 2pns reso- 
nance series stretches between the 2s threshold and the 
2p threshold, a range of about 4 eV. Our attention is di- 
rected at the time scale and energy scale of this series. 
The initial continuum distribution is set to be centered 
at 2.3 eV above the ionization threshold, with bandwidth 
of 1.825 eV, or pump duration of 1 fs. The bandwidth, 
covering the resonances of interest, is limited to the en- 
ergy range above the 2s threshold to avoid unnecessary 
complications. Each resonance term in Eq. [21] is multi- 
plied by a Gaussian window whose width is 2r to narrow 
down the interaction V to some finite range in the cal- 
culation. The resultant resonance profiles arc displayed 
in Fig. |6l along with the experimental photoabsorption 
cross section values in Ref. [22| > plotted by the light gray 
curve and normalized to fit in the vertical range of our 
calculation. The 2pns resonance series is seen in the ex- 
periment from low to high energy, starting with 2p3s at 
about 1.5 eV. For the calculation, the figure shows how 
the profiles evolve in time and are finalized, where the 
lower (wider) resonances build up quickly and the higher 
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FIG. 5: The resonance profile |d B /(f)| 2 for the 2p4s autoion- 
ization in Be, at (a) t — 0, 1.5 fs, 3 fs, and 4.5 fs, and at (b) 
t = 10 fs, 20 fs, and 40 fs. 



(narrower) resonances slowly. At short times, the contin- 
uum must be viewed as a whole distribution which gives 
no information on individual resonances; only at large 
times, e.g. t = 160 fs, the profiles assigned to individual 
resonances can be seen. Nonetheless, our calculation and 
the experiment in Ref. (22j should not be compared with 
each other directly. In the present work, a short pump 
pulse covers a wide energy range; the system propagates 
coherently with time. On the other hand, the experiment 
in Ref. [22j was done with many independent incoherent 
measurements, where each measurement was made for a 
single energy point using a long pulse of at least of hun- 
dreds of picoseconds, and gives back a single quantity for 
the cross section. 

We caution the readers that the time-dependent res- 
onance profile shown in Fig. [5] is from theoretical cal- 
culations only, i.e., a plot of |c?£;(i)| 2 vs E for different 
time delays after the pump pulse. As emphasized earlier, 
|d£;(£)| 2 cannot be measured directly in the laboratory 
since any physical measurement will take at least many 
picoseconds which cannot distinguish time evolution that 
happens on the femtosecond scale. Secondly, the wave 
functions associated with the coefficients de(t) in Eg. IT51 
are not eigenstates of the atom such that | cZ ^ (*) | 2 cannot 



be directly interpreted as measurable probabilities. 
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FIG. 6: The calculated resonance profiles for the 2pns ( P°) 
resonances at t = 0, 4 fs, 8 fs, and 160 fs in Be, and the 
experimental phot oabsorpt ion cross section in Be between the 
2s and the 2p thresholds [S^]. The experimental values are 
normalized to the vertical range of the calculation. 



IV. METHOD OF PROBING THE 
TIME-DEPENDENT AUTOIONIZATION 
PROCESSES 

After we have predicted how the wave packet in the 
autoionization process evolves with time, an experimen- 
tal scheme where such predictions can be confirmed is 
essential. Using SAP-pump and IR-probe, the lifetime of 
Auger decay of an inner-shell hole Q or of a Fano reso- 
nance in helium [l2| have been " deduced" . The streaking 
of an autoionizing electron in the IR laser field has been 
calculated by Wickenhauser et al ■ The streaked elec- 
tron spectra are too complicated to allow the retrieval of 
the predicted intermediate profiles. Taking beryllium as 
example, here we propose in the following a way to mea- 
sure the time dependence of Fano profiles. 

The probe scheme is sketched in Fig. [7] After the 
pump, we define the starting time for the autoioniza- 
tion as t = such that the pump field is finished for 
t > 0. A probe pulse is shined onto the system with 
its peak at time delay r. The probe has the duration 
of 1.5 fs, or bandwidth of 1.216 eV, and carrier energy 
uj of 40 eV. This probe pulse ionizes the 2s electron in 
the 2sEp continuum (binding energy of 2s is 18.21 eV). 
After the probe pulse is over, the excited wave packet of 
Eq. [30] changes to 

!*«*(*)} =XX(*)|2pns)+ f d E (t)\2sEp)dE 




d E , E {t)\E'pEp)dE'dE. (32) 
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We have assumed that the 2pns bound states are not 
changed by the probe pulse, and only a small part of 
the 2sEp continuum states are ionized (thus the ampli- 
tude now given by d E (t) is slightly different from d E (t)) 
to generate new E'pEp states. After the probe pulse, 
both d n {t) and d E (t) continue to change with time as 
the bound states decay. The newly created E'pEp states 
are eigenstates of the atomic (field-free) Hamiltonian and 
are stationary states, and they can thus be probed by a 
laboratory detector. Using 40 eV probe pulses, the en- 
ergy E' will be centered around 21.8 cV, well-separated 
from the electron energy of 2.5 eV due to autoionization. 
Clearly the probability \dE>E(t)\ 2 should be independent 
of time and is proportional to |d_E(r)| 2 . To determine 
\d E {T)\ 2 one can measure the energies of the two electrons 
in coincidence or by measuring the low-energy electrons 
in coincidence with Be 2+ ions. Either way lower energy 
group duplicates the original continuum profile at the 
moment of the probe. By changing the delay time of the 
probe a copy of |d£;(r)| 2 can be determined. Due to the 
finite duration of the probe pulse, the measured copy of 
\d E {r)\ 2 will be somewhat smoothed out. 



perturbation theory by 
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photocopy 
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g e 2+ no change [j e 2+ 
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FIG. 7: A conceptual scheme of probing the time-dependent 
resonance profiles. At t = the pump pulse is over and au- 
toionization starts while the system evolves freely. To find 
out the the evolving electron wave packet at time r, a probe 
pulse is used to ionize the 2s electron in the 2sEp continuum 
and create the E'pEp states. After the probe is over, the new 
wave packet consists of a part of the old wave packet which 
continues to evolve and autoionize (the upper row) and the 
newly generated part made of E'pEp states (the lower row). 
The latter does not change with time. By picking the respon- 
sible energy range and collecting the double ionized signals, 
the resonance profile in the autoionization at the moment r 
is captured. 



If the probe is of perturbative strength, we can calcu- 
late the predicted two-electron spectra quantitatively. In 
the presence of the external probe field i^W, the coeffi- 
cient for E'pEp is given by the first-order time-dependent 



d E > E = ( E'pEp 



U 



To To 
2 ' 2 



^ex ~ 



T 



(33) 



111 



111 



E'pEp |e- 4/ M^-*V 2 (i)| dt 

(34) 

F 2 (t)d E >(T + t)dt (35) 



where V 2 {t) is the dipole interaction potential for the 
probe pulse. The interaction is turned on only for 
—To/2 < t < To/2. Note that t is now measured with 
respect to the center of the probe, and the total ex- 
cited wave function \^ ex (—To/2)) at t = —To/2 has al- 
ready gone through the autoionization process for a du- 
ration of t — To/2. In Eq. I33ll3"4"l the time-propagator 
cxp [— iHo(To/2 + t)] in U operates on ^ ex (—To/2) and 
gives \& e:c (t) in return. In the next step in Eq. 1351 since 
the bound states 2pns do not participate in the probe 
process, only the continuum part is written down. The 
dipole matrix element p, = (E'p\z\2s) is assumed to be 
a constant. The polarization of the probe field is in the 
z-direction, i.e. F 2 (t) = zF 2 (t). The retrieved profile 
as a function of E is defined as J \d E > E \ 2 dE' . Fig. [5] 
shows very good agreement between the retrieved profile 
at different r's and the original continuum profile at the 
associated time r. 
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FIG. 8: The time-dependent continuum profile compared with 
the probed profile in the Be case, (a) The calculation as spec- 
ified in Fig. [6] (b) The retrieved profile using a probe beam 
with to = 40 eV and the duration (FWHM) of 1.5 fs. The 
signal values have been scaled to match the original profile at 
t — 160 fs. The agreement between (a) and (b) is good. 

While the calculation given above is in first-order per- 
turbation, the scheme can be applied to intense probe 
as well as long as the probe field dominantly removes 
the 2s electron and has little effect on other components 
in the system. The Ep electron is not affected by the 
probe, and the continuum profile of 2sEp at r will still 
be transferred to E'pEp plainly (with additional peaks 
possibly from two-photon ionization). The advantage of 
a stronger probe is that it increases the signal rate. This 
may be an important factor to realize this pump-probe 
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scheme since the population of double excited states by 
the pump can be already weak. In the present example, 
the photon energy is 11.62 eV for the pump and 40 eV 
for the probe, presuming that the short UV sources are 
provided for both the pump and the probe. 

V. CONCLUSIONS 

In this paper we study the temporal behavior of the 
atomic autoionization process. Following the original pa- 
per of Fano, we treat autoionization as the decay of a dis- 
crete state into a degenerate continuum. Fano's theory 
characterizes such a resonance with a shape parameter 
q and a width T, in addition to the resonance energy 
position. For many decades, Fano resonances are mea- 
sured and characterized in energy domain experiments. 
With the advent of XUV and soft X-rays with pulse du- 
rations of attoseconds to a few femtoseconds, we address 
the question whether Fano resonances can be probed di- 
rectly in the time domain within its typical lifetime. In 
particular, can we talk about the time-dependent Fano 
profiles before the discrete state is fully decayed and can 
such profiles be measured experimentally? 

To describe time-dependent Fano profiles, it is inconve- 
nient to express the total time-dependent wave function 
in terms of the eigenstates of the system. Instead, the 
configuration basis states, including the bound and the 
continuum, used in Fano's original theory are much more 
useful (or more physical) for analyzing the autoionization 
process, and we are led to define the time-dependent Fano 



profiles by referring to the wave packets expanded in the 
continuum configuration basis states. To probe the tem- 
poral Fano profiles with good time resolution, we suggest 
that the decay process be perturbed by a short soft-X-ray 
pulse which interrupts the autoionization at the time of 
the probe. The probe in the meantime creates a station- 
ary wave packet that can be measured in the laboratory. 
We further consider many Fano resonances together if 
they fall within the bandwidth of the pump pulse. 

The pump-probe scheme proposed here, requiring 
high-energy short pulses for both the pump and the 
probe, is not yet available today. On the other hand, 
attosecond technology is undergoing rapid growth, the 
proposed scheme of measuring autoionization dynamics 
may become doable within a few years. At the same time, 
theory for attosecond electron dynamics is still a bar- 
ren field today. General approaches for understanding 
electron dynamics for many-electron systems are badly 
needed at this time. Not all the probe pulses have sim- 
ilar effectiveness in retrieving the time-dependent wave 
packets. Without the general theoretical framework, it 
would be difficult to retrieve electron dynamics from the 
signals measured by the probe pulse. 
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